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Abstract 

The statistical properties of dynamically triangulated manifolds (DT mfds) in 
terms of the geodesic distance have been studied numerically. The string suscep- 
tibility exponents for the boundary surfaces in three-dimensional DT mfds were 
measured numerically. For spherical boundary surfaces, we obtained a result con- 
sistent with the case of a two-dimensional spherical DT surface described by the 
matrix model. This gives a correct method to reconstruct two-dimensional random 
surfaces from three-dimensional DT mfds. Furthermore, a scaling property of the 
volume distribution of minimum neck baby universes was investigated numerically 
in the case of three and four dimensions, and we obtain a common scaling structure 
near to the critical points belonging to the strong coupling phase in both dimen- 
sions. We have evidence for the existence of a common fractal structure in three- 
and four-dimensional simplicial quantum gravity. 
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1 Introduction 



There has been, over the last few years, remarkable progress in the quantum theory of 
two-dimensional gravity. Two distinct analytic approaches for quantizing two-dimensional 
gravity have been established. These are recognized as a discretized theory [0] and a con- 
tinuous [@] theory. The discretized approach, implemented by the matrix model technique, 
exhibits behavior found in the continuous approach, given by Liouville field theory, in a 
continuum limit, for example, the string susceptibility exponent and Green's function. 
Strong evidence thus seems to exist for the equivalence between the two theories in two 
dimensions. A numerical method based on the matrix model, such as that of dynamical 
triangulation [Q], has drawn much attention as an alternative approach to studying non- 
perturbative effects, also being capable of handling those cases where analytical theories 
cannot yet produce meaningful results. In the dynamical-triangulation method, calcu- 
lations of the partition function are performed by replacing the path integral over the 
metric to a sum over possible triangulations. Over the past few years a considerable 
number of numerical studies have been made on three- and four-dimensional simplicial 
quantum gravity [|J, and recent results obtained by dynamical triangulation in three and 
four dimensions suggest the existences of scaling behaviors near to the critical point |5],|J 
in terms of the geodesic distance 0. When the coupling strength (k) becomes close to the 
critical point from the strong-coupling side (k < k c ), which corresponds to the crumple 
phase, the transition is smooth. On the other hand, when k becomes close to the critical 
point from the weak-coupling side (k > k°), which corresponds to the branched-polymer 
phase, the transition is very rapid. In fact, in ref. H it is reported that near to the critical 
point belonging to the strong-coupling phase (k — > k c from k < k c ) the scaling behavior 
of the mother boundary surfaces exists at appropriate geodesic distances, and also that 
no mother boundary surface exists in the weak-coupling phase (i.e., branched polymer 
phase). It seems reasonable to suppose that the model makes sense as long as we are 
close to the critical point from the strong-coupling phase. More noteworthy is that there 
are several pieces of evidence that the phase transition of three- and four-dimensional sim- 
plicial gravity is first order |8|,|9|]. We actually observe double-peak histogram structures 
which are signals of a first-order phase transition for appropriate large lattice sizes. There- 
fore, we carefully chose one of the peaks belonging to the strong phase as an ensemble for 
our simulations. 

We must thus look more carefully into these boundary mfds. In the last few years, 
several articles have been devoted to the study of boundary mfds. In two dimensions, it is 



revealed by ref. |10[ that the dynamics of the string world sheet (random surfaces) can be 
described by the time (i.e., geodesic distances) evolution of boundary loops. Ref. JXTJ] is 



also based on the idea that the functional integral for three- and four-dimensional quantum 
gravity can be represented as a superposition of a less-complicated theory of random- 
surface summation over a two-dimensional surface 0. It is precisely on such grounds 
that we claim that the higher dimensional complicated theory of quantum gravity can be 
reduced to the lower dimensional quantum gravity. 

This paper is organized as follows. In Sec. 2 we briefly review the model of the three- 
dimensional dynamical triangulation. In Sec. 3 we report our measurements for the string 



*In this case a definition of a time direction is different from ours. 
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susceptibility exponent (7^) of the mother boundary surfaces in three-dimensional DT 
mfds. In Sec. 4 we report on our measurements for the string susceptibility exponent (7^) 
of the three- and four-dimensional DT mfds. We summarize our results and discuss some 
future problems in the final section. 

2 Model 

We start with the Euclidean Einstein-Hilbert action in D(= 3,4) dimensions, 

S EH = Jd D x,/g^A-^Rj, (1) 

where A is the cosmological constant and G is Newton's constant of gravity. We use the 
lattice action of the D-dimensional model with the S D topology, corresponding to the 
above action, as follows: 

S(k , k d ) = -k N + K D N D 

= "§iVo+(A , -^(27r-6cos- 1 (i)))iV D , (2) 

where Ni denotes the total number of i-simplexes, and A' = cA; c is the unit volume of a 
D-simplexes, and cos _1 (-^) is the angle between two tetrahedra. The coupling («o)S is pro- 
portional to the inverse of the bare Newton's constant, and the coupling (kd) corresponds 
to a lattice cosmological. 

For the dynamical triangulation model of D-dimensional quantum gravity we consider 
a partition function of the form 

Z(ko,k d )= £ e -s(^). (3) 

T(S D ) 

We sum over all simplicial triangulations (T(S D )) on a D-dimensional sphere. In practice, 
we must add a small correction term (AS) to the lattice action in order to suppress any 
volume fluctuations. The correction term is denoted by 

AS = 5(N D - A^ ar S et) ) 2 , (4) 
where _/v^ ar S e ^) j s fa e target volume of D-simplexes, and we use S = 0.0005 in all runs. 

3 Boundary structure of three-dimensional DT mfds 

We now define the intrinsic geometry using the concept of a geodesic distance as a min- 
imum length (i.e., minimum step) in the dual lattice between two 3-simplexes in a three 
dimensional DT mfd. Suppose a three-dimensional ball (3-ball) which is covered within d 
steps from a reference 3-simplex in the three-dimensional mfd with S 3 topology. Naively, 

tin the case of four dimensions, we use K2 in all runs instead of k with the conventions. 
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Figure 1: Schematic picture of boundary surfaces (Si,E2 and S3) at distance d and a 
surface (S' 2 ) at distance d! in 3D Euclidean space M with S 3 topology. The mother 
boundary surface is defined as a surface (S 2 ) with the largest tip volume (V2), and the 
other surfaces are defined as the baby boundary. 



the 3-ball has a boundary with spherical topology (S 2 ). However, because of the branch- 
ing of DT mfds, the boundary is not always simply-connected, and there usually appear 
many boundaries which consist of closed and orientable two-dimensional surfaces with any 
topology and nontrivial structures, such as links or knots i. Below, we consider only three- 
dimensional DT mfds with S 3 topology. We can show a sketch of typical configuration 
in Fig.|I]. M denotes a three-dimensional DT mfds with S 3 topology and Si, £2 and S3 
denote the boundaries which are closed and orientable triangulated surfaces at a distance 
d from an origin (x). In a previous paper [T2|] we showed that the coordination number 



(q) distributions of the spherical mother boundary surfaces of three-dimensional DT mfds 
are consistent to the theoretical prediction for the two-dimensional random surfaces in 
the large-g region. 

Here, we concentrate on the boundary surfaces in three dimensions. The ensemble of 
these boundary surfaces consists of various volume surfaces. The boundary surfaces are 
divided into two classes: one is a baby boundary and the other is a mother one. Here, 
we give a precise definition of "baby" and "mother" universes (see Fig.|I|). The mother 
universe is defined as a boundary surface (S 2 ) with the largest volume (V2), and the 
other surfaces are defined as baby universes. The boundaries with small sizes (~ 0(1)) 
are called a baby one. The baby boundaries originate from the small fluctuations of the 
three-dimensional Euclidean spaces. We thus think that these surfaces are non-universal 
objects. In ref. || it is reported that the surface-area-distributions (SAD) of the mother 
universe show the non-trivial scaling behavior near to the critical point. If this scaling 



tin the strict sense links or knots are constructed by loops. In our case the loop is a fat loop. 
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x = B(1-B/A) 



Figure 2: Volume distribution of minbu in terms of the spherical boundary surfaces in 
3D DT mfds at distance 10. Total volume of 3D DT mfds is 16K. The lower limits of the 
volume of the spherical boundary surfaces is 400. 

behavior turns out to be correct in the strict sense of the limit iV — > oo, we can take 
a continuum limit of three-dimensional DT mfds as well as two-dimensional DT mfds. 
Therefore, we focus on the mother boundary near to the critical point in the following. 

The following question now arises: how to measure -y s t of these boundary surfaces with 
different sizes. Suppose that a spherical boundary surface with area A is obtained; we 



can then use the standard minbu algorithm [13] in order to measure jst- The distribution 



function for the minbu with size B can be written as 

3(A-B + l)(g + \)Z\A -B + l]Z[B + 1] 
n A {B) = ^ . (5) 

When we substitute the asymptotic form of the partition function, Z[A] ~ j4 7s '~ 3 e MA , for 
eq.(||) we obtain the normalized distribution function, 

n A (B) = c{(B + l)(l-— — )} , (6) 

where c denotes a normalization factor which depends on the area (A). The spherical 
mother boundaries are selected as grand-canonical ensembles. In the case that the topol- 
ogy of the boundary surface is not a sphere, but a handle-body, such as a torus, the 
standard minbu algorithm does not work welll. It is for this reason that we concentrate 
on the spherical boundary surface. 

There are other things to note in terms of the grand-canonical ensembles in question. 
Using the grand-canonical ensembles for the measurement j st of boundary surfaces, we 



§To put it another way, when the topologies of the baby boundary and the mother boundary are 
different, the minbu algorithm does not work well. 



4 




10 

x = B(1-B/A) 



Figure 3: Volume distribution of minbu in terms of the spherical boundary surfaces in 
3D DT mfds at distance 11. Total volume of 3D DT mfds is 32K. There are several lower 
limits of the volume of the spherical boundary surfaces: 200, 400 and 600. 



suffer from finite-size effects. Therefore, we should introduce a lower limit of the sizes 
of the boundary surfaces. Then, the grand-canonical ensemble does not contain surfaces 
whose sizes are less than the lower limit. We choose an appropriate geodesic distance 
corresponding to the total volume (A^). 

Fig.0 shows the distribution of the minbu of the spherical mother boundary surfaces at 
d= 10 with a total volume (A3) of 16K. Fig.|3] shows the distribution of the minbu of the 
spherical mother boundary surfaces at d= 11 with A3 = 32K. Whether those introduced 
lower limits are sufficient or not is open to question, and the average size of the grand- 
canonical ensembles is about 700 for lower limit of 600 in Fig.|3|. We examed separately 
the string susceptibility exponent in two-dimensional pure gravity with a total area of 
700, and actually obtained 7^ = —0.49(5) with the number of ensembles being 10 5 , even 
in such a small-size simulation. Fig.|] reveals that the scaling property of the minbu area 
distributions becomes more clear the bigger is the average size of the boundary. Thus, 
these distributions remain after the thermal limit A3 — * oo. It should be concluded, 
due to these scaling relations, that the string susceptibility exponent of the spherical 
mother boundary surfaces in three dimensions is completely consistent with the string 
susceptibility exponent of the random surface in two dimensions. Our numerical results 
given here agree with the discussions in refs. |I1,[TH. 



On the other hand, when we apply the method in the previous section to higher 
dimensions, we are confronted by a difficulty. We cannot use Euler's character in order 
to distinguish the topologies of three-dimensional DT mfds. That is to say, to chose 
a spherical three-dimensional boundary mfd in four- dimensional mfd is not easy at all. 
Thus, it is difficult to determine j s t of the boundary mfds in four dimensions by means of 
the minbu technique because of the reason mentioned above. It is too involved a subject 
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Figure 4: Minbu volume distribution ua{B) of the three-dimensional DT mfds with 
N3 = 16K with log-log scales. Coupling constants kq = 4.09 and K3 = 2.20 for N3 = 16K. 
Dotted line is drown using the least squares with the data points of 800 configurations 
with a range of (x=)40 ~ 200. 

to be treated here in detail. 

Thus in the following section we consider the scaling structures of the minbu distribu- 
tions in three- and four- dimensional simplicial quantum gravity. 



4 Fractal structures of three- and four- dimensional 
DT mfds 

In the introduction we report on the possibility that this model in three- and four- 
dimensional simplicial quantum gravity makes sense near to the critical point belonging 
to the strong-coupling phase. Therefore, we investigated the scaling structures of the 
minbu distributions near to the critical point belonging to the strong coupling phase nu- 
merically. Fig.[| shows the distribution of the minbu of three-dimensional DT mfds with 
a total volume (N3) of 16K. We tune kq = 4.09 and K3 = 2.20. The three branches 
shown in Fig.[| reflect the symmetric factor of three-dimensional DT mfds. The influence 
of such a symmetric factor has not been observed in two-dimensional DT surfaces. Fig.[5] 
shows the distribution of minbu of four-dimensional DT mfds with a total volume (N4) 
of 6AK. We tune k% = 1.275 and K4 = 1.353. We also observe two branches which 
reflect the symmetric factor of four-dimensional DT mfds. In three and four dimensions 
we have a similar scaling behaviour of the minbu volume. From these scaling date which 
are shown in Figs.|, || we can extract the string susceptibility exponents: j st = —0.05(6) 
with N 3 = 1QK and j st = 0.26(5) with iV 4 = UK pf. 
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Figure 5: Minbu volume distribution ua{B) of the four-dimensional DT mfds with 
N± = 64K with log-log scales. Coupling constants K2 = 1.275 and K4 = 1.353 for 
N4 = 64K. Dotted line is drown using the least squares with the data points of 500 
configurations with a range of (x=)40 ~ 1000. 

4.1 Continuum limit in four-dimensional pure simplicial quan- 
tum gravity 

It is generally agreed that the phase transition of simplicial quantum gravity in both three 
and four dimensions is first order |8],|9[] . Why is the scaling behaviour observed near to the 
critical point? Let us devote a little more space to discussing this question. In ref. |L7| a 
quantum theory of four- dimensional gravity which is the analog of the Liouville theory of 
two-dimensional quantum gravity has been argued by considering the effective action for 
the conformal factor of the metric induced by the four-dimensional trace anomaly. The 
authors in the reference have also argued for a model including the volume (cosmological 
constant) and Einstein terms. The string susceptibility exponent (7^) obtained by ref. [[[7] 
in four dimensions has a similar behaviour to the j st obtained by the Liouville theory in 
two dimensions: 



% t (Q 2 ) = -^(1-7^ + ^-7^) ( 4Dim ) (7) 

and 



7«t(cj= |2 



l(l- C+v /(l-c)(25-c)) (2 Dim). (8) 



Here, Q 2 which plays the role analogous to matter central charge in two dimensions is 
the coefficient of the Gauss-Bonnet term in the trace anomaly []T7|]. If Q 2 > 8 (i.e., 



lst{Q ) < 0, the cosmological and the Einstein terms become irrelevant in the infinite 
volume limit, and the spike configurations are suppressed [II]]. Smooth configuration (i.e. 



smooth phase) is expected to appear. Q 2 = 8 (7^ = 0) behaves in many respects like the 
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Figure 6: Phase diagram expected from numerical simulations for pure 4D simplicial 
gravity, df indicates the fractal dimension obtained in 4D. 



c = 1 case in two-dimensional quantum gravity. 

Our numerical results mentioned in the previous subsection, as shown in Fig|| suggest 
that there are no signals for the smooth phase in which we can take the continuum limit. 
It seems reasonable to suppose that the scaling relations which we obtained near to the 
critical point are pseudo-scaling relations. 

5 Summary and Discussion 

We investigated the statistical properties of boundary surfaces of three-dimensional DT 
mfds with S 3 topology near to the critical point. We have found positive evidence that the 
spherical mother boundary surfaces in three dimensions are equivalent to two-dimensional 
spherical random surfaces described by the matrix model. We thus have a conjecture: the 
mother boundary surfaces in three dimensions with any handles are statistically equiv- 
alent to two-dimensional random surfaces described by the matrix model. Therefore, if 
these boundary surfaces in three dimensions can be recognized as random surfaces, three- 
dimensional DT mfds can be reconstructed by the direct products of X (two-dimensional 
DT surfaces) and d (geodesic distance)!. Unfortunately, we cannot use the Euler's char- 
acter in order to distinguish the topologies of three-dimensional DT mfds. Thus, it is 
difficult to determine j st of the boundary mfds in four dimension by means of the minbu 
technique. There is room for further investigation. 

Furthermore, we investigated the volume distributions of the minbu of three- and four- 
dimensional DT mfd near to the critical point. Our numerical results of the minbu volume 
distribution in three and four dimensions show a tendency which is also similar to two- 
dimensional case. If we assume that the partition function in three and four dimensions 

^This strategy can straightforward be applied to the four-dimensional case if the equivalence between 
the boundary DT mfds in four dimensions and the three-dimensional DT mfds is firmly established. 
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can be written in the same asymptotic form as the two-dimensional case, we obtain the 
string susceptibility exponents: ^ st = —0.05(6) with N 3 = 16K and j st = 0.26(5) with 
A^ 4 = 64K near to the critical point. 
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